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Abstract 

We consider the possibility that the gravitino might be the lightest supersymmetric par- 
ticle (LSP) in the constrained minimal extension of the Standard Model (CMSSM). In this 
case, the next-to-lightest supersymmetric particle (NSP) would be unstable, with an abun- 
dance constrained by the concordance between the observed light-element abundances and 
those calculated on the basis of the baryon-to-entropy ratio determined using CMB data. 
We modify and extend previous CMSSM relic neutralino calculations to evaluate the NSP 
density, also in the case that the NSP is the lighter stau, and show that the constraint from 
late NSP decays is respected only in a limited region of the CMSSM parameter space. In 
this region, gravitinos might constitute the dark matter. 
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1 Introduction 



If R parity is conserved, the lightest supersymmetric particle (LSP) is stable, and a possible 
candidate for the cold dark matter postulated by astrophysicists and cosmologists [1] . Most 
analyses of such supersymmetric dark matter have assumed that the LSP is the partner of 
some combination of Standard Model particles, such as the lightest neutralino Xi with an 
abundance calculated from the freeze-out of annihilation processes in a thermal initial state. 
However, another generic possibility is that the LSP is the gravitino G [2] - [8], whose relic 
abundance would get contributions from the decays of the next-to-lightest supersymmetric 
particle (NSP) and possibly other mechanisms. 

As we discuss in more detail below, the lifetime of the NSP is typically such that it decays 
between Big-Bang nucleosynthesis (BBN) and the 're-' combination process when the cosmic 
microwave background (CMB) was released from matter. Since NSP decays release entropy 
during this epoch, they are constrained by the concordance of the observed light-element 
abundances with BBN calculations assuming the baryon-to-entropy ratio inferred from CMB 
observations. For a typical lifetime t^sp = 10^ s, the observed ^Li abundance implies [9] 

< 5 X io~" fl5i^) (1) 

V niNSP / 

before NSP decay, with the D/H ('^He) abundance providing a constraint which is weaker 
by a factor of about 10 (20). Assuming a baryon-to-entropy ratio r] = nB/n^ = 6.0 x 10^^*^, 
in agreement with the WMAP result i] = 6.l1lo[2 ^ 10^^° [10], implies the constraint 
nNSp/^B < 10^''(100 GeV/niNSp) before the onset of NSP decay. To assess the power of 
this constraint, we re-express it in terms of fl%gph'^, the relic density that the NSP would 
have today, if it had not decayed: 

fl%sph'^ < 10-2 Qsh'^ ^ 2 X 10-^ (2) 

where Qsh'^ ~ 2 x 10^^ is the present-day baryon density. However, the requirement (j2I) 
would be relaxed for a shorter-lived NSP [7], as we discuss later. 

In contrast, assuming that the lightest neutralino x is the LSP, there have been many 
calculations of f^^^^ constrained minimal supersymmetric extension of the Standard 

Model (CMSSM), in which the GUT-scale input gaugino masses mi/2 and scalar masses mg 
are each assumed to be universal [11-14]. These calculations find generic strips of CMSSM 
parameter space in which 

n^h^ ~ 5 X Qeh^ ~ 0.1 (3) 



1 



This is similar to the range of the cold dark matter density QcDMh"^ favoured by astro- 
physicists and cosmologists, which is one reason why neutralino dark matter has been quite 
popular. 

In this paper, we assume no a priori relation between m3/2 and the soft supersymmetry- 
breaking masses mi/2 and mo of the spartners of Standard Model particles in the CMSSM. 
This is possible in, e.g., the framework of = 1 supergravity with a non-minimal Kahler 
potential [15]. In such a framework, the LSP might well be the gravitino G. In this case, the 
NSP would likely be the lightest supersymmetric partner of some combination of Standard 
Model particles, such as the lightest neutralino x o^' the lighter stau fi. Particularly in the 
X NSP case, one might expect fl%gph'^ to be near the range Comparing this with the 
condition Q necessary for gravitino dark matter, we see that, if t^vsp = 10^ s, gravitino dark 
matter could be possible only in rather different regions of the CMSSM parameter space, 
where the NSP density is very suppressed compared with the usual x density. Moreover, 
in this case, NSP decays alone could not provide enough gravitinos, since they could only 
yield ^^3/2/^^ < ^atsp^^; there would need to be some supplementary mechanism for 
producing gravitinos, if they were to provide all the cold dark matter. For example, gravitino 
production during reheating after inflation could produce a sufficient abundance of gravitinos 
if the reheat temperature is relatively large, ~ 0(10^'') GeV [7]. 

The ffist step in our exploration of the gravitino dark matter possibility is to calculate 
n^5p/i^ throughout the {mi/2,Tno) planes for different choices of tan/? and the sign of fi in 
the CMSSM, assuming that the trilinear soft supersymmetry-breaking parameter Aq = 0. 
In the regions where < m^^, this is essentially equivalent to the usual neutralino dark 
matter density calculation. However, as we discuss below, this calculation must be adapted 
in the region where m^^ < m^. Moreover, one must take into account the possibility of 
a cosmological fi asymmetry, in which case the relic fi density would be larger than that 
given by the standard freeze-out calculation. We next compute the NSP lifetime and use 
the detailed constraints from the abundances of the light elements as computed in (0) for 
fixed ?7 = 6 X 10^^°. This allows us to delineate the regions of the CMSSM (mi/2, ^o) planes 
where gravitino dark matter appears possible. We find limited regions of the (mi/2,'^0) 
planes that are allowed. In these regions, the density of relic gravitinos due to NSP decay 
is typically less than the range favoured by astrophysics and cosmology. As noted above, 
supplementary mechanisms for gravitino production, such as thermal production in the early 
Universe, might then enable gravitinos to constitute the cold dark matter. 
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2 NSP Density Calculations 



In the framework of the CMSSM with a Ught gravitino discussed here, the candidates for 
the NSP are the hghtest partners of Standard Model particles. In generic regions of CMSSM 
parameter space, these are the lightest neutralino x ^-^^d the lighter stau fi ^. In regions 
where x is the NSP, the calculation of the NSP density ^%sp^^ identical with that of 
VLisph"^ in the CMSSM with a heavier gravitino, and we can recycle standard results. 

Extending these calculations of ^%sp^^ ^o regions where the fi is the NSP requires some 
modifications. Whereas the Majorana x is its own antiparticle, one must distinguish between 
the fi and its antiparticle f*, and calculate the sum of their relic densities. This requires 
a careful accounting of the statistical factors in all relevant annihilation and coannihilation 
processes. We have also made a careful treatment of the regions where there is rapid fi — rj* 
annihilation via Higgs poles, and a non-relativistic expansion in powers of the NSP velocity 
is inadequate. Here our treatment follows that of the neutralino LSP case in [11, 13, 17]. 

It is important to note that one would, in general, expect a net fi asymmetry rjfj^ = Xtjb, 
where A ~ 0(1). This would be the expectation, for example, in leptogenesis scenarios, 
and would also appear in other baryogenesis scenarios, as a result of electroweak sphalerons. 
However, in the context of the MSSM, there exist fifi — > rr annihilation processes which 
would bleed away any existing lepton asymmetry stored in the f sleptons, and the final relic 
density is given by the calculation described above. 



Using the standard = 1 supergravity Lagrangian [18,19], one can calculate the rates for 
the various decay channels of candidate NSPs to gravitinos. 

The dominant decay of a x NSP would be into a gravitino and a photon, for which we 
calculate the width 



where C^^ = {Oi^cos 6w + 02x^^^(^w) and O is the neutralino diagonalization matrix, 

M.nO = A4'j^"'^ . Note that in this and the following equations Mp = 1 /^/StiGn- 

^The lighter stop could also be the NSP if the trilinear coupling Aq is large [16], but here we fix = 
for simplicity. 



3 NSP Decays 
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A X NSP may also decay into a gravitino and a Z boson, for which we calculate the rate 

1 

^X^GZ = 772" 2^-^("^X''^3/2,Mz) 



X 




/I ml/A Ml 



where C^z = {.~Oix sin Ow + ^2;^ cos6'vy), and we use the auxihary functions 

2\ / / , ^ \ 2~ 



J='{m^, 1713/2, Mz] 



^ my2 + Mz \ \ / 7713/2 -M, 



771^ / I \ \ 771^ 



1/2 



(6) 



^K,m3/2,M,) = 1 + ^ (^-4 + ^j+^- — . (7) 

Note that in the hmit Mz — > we obtain T^^q^ '^x-^G-y ^"7 replacing C-^z with C^^. 
Decays of a x NSP into a gravitino and a Higgs boson are also possible, with a rate 

1 C^^j 771^ 

771^ ) 77llj \\77l^ J 77ll 

where C^h = {O^^ cos a — O^^ sin a). 

Analogously, for the heavy Higgs boson H we get T^^g/i ~^ ^x^gh replacing C^h with 
C^H = {O^^ sin a + O^^ cos a), and 77ih with m//. The corresponding formula for x ^ G + A, 
where A is the CP-odd Higgs boson in the MSSM, is also given by (jH)), but with 77ih replaced 
by 7nA, 7713/2 -7713/2 and C^h ^ C^a = (O^^ cos/5 + 03^ sin/5). 

Finally, the dominant decay of a f NSP would be into a gravitino and a r, with the rate: 



771^ ( "^3/2 V 



ip /;t3/2 

where we have neglected the O {m^ / tu^^) terms. 

4 Effects of Gravitino Decay Products on Light-Element 
Abundances 

The effects of electromagnetic shower development between Big-Bang Nucleosynthesis (BBN) 
and 're-'combination have been well studied, most recently in [9], where the simplest case of 
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X ^ G + 7 decays were considered. The late injection of electromagnetic energy can wreak 
havoc on the abundances of the light elements. Energetic photons may destroy deuterium, 
destroy "^He (which may lead to excess production of D/H), destroy ''Li, and/or overproduce 
^Li. The concordance between BBN calculations and the observed abundances of these ele- 
ments can be used to derive a limit on the density of any decaying particle. In general, this 
limit will depend on both the baryon asymmetry rjB, which controls the BBN predictions, 
and on the life-time of the decaying particle tx- For a fixed value r^^ = 6 x 10~^°, as sug- 
gested by CMB observations, the bounds derived from Fig. 8(a) in [9] may be parameterized 
approximately as 

?/< 0.13 - 2.85 x + 3.16, (10) 

where y = log(Cx/GeV) = log(mx^js:/'"-7/GeV) and x = log(rx/s), for the electromagnetic 
decays of particles X with hfetimes 10^^ s > t^sp ^ 10'^ s. In our subsequent analysis we 
use the actual data corresponding to the limit in [9] in order to delineate the allowed regions 
of the (mi/2, mo) planes, but (fTUj) may help the reader understand qualitatively our results. 

The other NSP decay modes listed above inject electrons, muons and hadrons into the 
primordial medium, as well as photons. Electromagnetic showers develop similarly, whether 
they are initiated by electrons or photons, so we can apply the analysis of [9] directly also to 
electrons. Bottom, charm and r particles decay before they interact with the cosmological 
medium, so new issues are raised only by the interactions of muons, pions and strange 
particles. In fact, if the NSP lifetime exceeds about 10^ s, these also decay before interacting, 
and the problem reduces to the purely electromagnetic case. In the case of a shorter-lived 
NSP, we would need to consider also hadronic interactions with the cosmological medium [20], 
which would strengthen the limits on gravitino dark matter that we derive below on the 
basis of electromagnetic showers alone. In the following, we do not consider regions of the 
(mi/2,mo) planes where t^sp < 10^ s. 

It is sufficient for our purposes to treat the decays of /i, vr and K as if their energies were 
equipartitioned among their decay products. In this approximation, we estimate that the 
fractions of particle energies appearing in electromagnetic showers are vr" : 100%, fi : 1/3, 
TT^ : 1/4, : 0.3, : 0.5. Using the measured decay branching ratios of the r, we then 
estimate that ~ 0.3 of its energy also appears in electromagnetic showers. In the case of 
generic hadronic showers from Z or Higgs decay, we estimate that ~ 0.6 of the energy is 
electromagnetic, due mainly to 7r° and vr^ production. 

Our procedure is then as follows. First, on the basis of a freeze-out calculation, we 
calculate the NSP relic density = 3.9 x 10^ GeV"^ (x- Next, we use the calculated 

life-time tx to compute the ratio of the relic density to the limiting value, Cx^^^ provided by 
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the analysis of [9] , taking into account the electromagnetic energy decay fractions estimated 
above. Finally, we require 



Cx 



5 Results 



As compared to the case of CMSSM dark matter usually discussed, in the case of gravitino 
dark matter one must treat m3/2 as an additional free parameter, unrelated a priori to ttlq 
and mi/2. We incorporate the LEP constraint on rrih in the same way as in [13] ^, and it 
appears as a nearly vertical (red) dot-dashed line in each of the following figures. Regions 
excluded by measurements of 6 — > 37 are shaded dark (green). For reference, the figures 
also display the strips of the (mi/2,'^0) planes where 0.094 < ^%sp < 0.129. This density 
is the same as Qisph"^ in a standard CMSSM analysis with a heavy gravitino, extended to 
include the unphysical case where the fi is the LSP. We note the familiar 'bulk' regions and 
coannihilation 'tails', as well as rapid-annihilation 'funnels' for large tan/? [17,21]. If these 
figures were extended to larger mo, there would also be 'focus-point' regions [22,23]. 

We now summarize our principal results, describing the interplay of these constraints 
with those associated specifically with gravitino dark matter, studying the {mi/2,fnQ) planes 
for three choices of tan/3 and the sign of fi: (1) tan/3 = 10, yU > 0, (2) tan/3 = 35, /i < 0, and 
(3) tan/3 = 50, yU > 0. In each case, we consider four possibilities for two fixed values 

10 GeV and 100 GeV, and two fixed ratios relative to mo: m3/2 = 0.2 mo and mo itself. If 
m3/2 ^ rrio, the G is typically not the LSP, and this role is played by the lightest neutralino 
X, as assumed in most analyses of the CMSSM. In each (mi/2,'n^o) plane, we display as a 
(purple) dashed line the limit where the density of relic gravitinos from NSP decay becomes 
equal to the highest cold dark matter density allowed by WMAP and other data at the 2-a 
level, namely ^^3/2/^^ < 0.129: only regions below and to the right of this contour are allowed 
in our analysis. 

Fig. [T] displays the (mi/2,'"^o) planes for tan/3 = 10 and yU > ^. Panel (a) displays the 
choice m3/2 = 10 GeV, in which case the LSP is the G throughout the displayed region of the 
(mi/2, fno) plane. Above and to the left of the (purple) dashed line, the relic density ^3/2/^^ of 
gravitinos yielded by NSP decay exceeds the 2-a upper limit on the cold dark matter density, 
0.129, imposed by WMAP and other cosmological data. This region is therefore excluded. 

^For simplicity, we do not show the LEP constraints on m^± and mg, which do not impinge on the regions 
of parameters aUowed by other constraints. 

■^The case tan/5 = 10 and < is very similar, with the exception that the b ^ constraint is more 
important. 
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In the regions below the (purple) dashed line, the relic G density might be increased so as 
to provide the required cold dark matter density if there were significant thermal gravitino 
production, in addition to that yielded by NSP decay. 

The light-element constraint on NSP decays is shown as the grey (khaki) solid line cor- 
responding to r = 1, where r is defined in Regions to the right and below this line 
are allowed by this constraint. Here, and in the remaining figures below, the region which 
satisfies the abundance constraint is labelled r < 1. There is a solid (black) line with 
mi/2 ~ 800 GeV which indicates where t^sp = 10"^ s. To the right of this line, t^sp < 10^ s, 
the case we do not consider here because additional constraints due to hadronic decays must 
be included, so this region is left blank ^. Here and in subsequent figures, the region that is 
allowed by all the constraints is shaded in light (yellow) color. 

We see that there is an extended strip between the grey (khaki) solid line and the solid 
(black) line. This strip is truncated above mo — 650 GeV, because the relic density of 
gravitinos from NSP decay becomes too large. This is true up to ~ 2900 GeV, where the 
relic density drops as we approach the focus-point region. Here a small allowed region opens 
up as the r = 1 curve bends towards lower values of mi/2- The allowed strip broadens in the 
low-mo region where m^^ < m^, below the dotted (red) line where m^ = m^^. In this region, 
gravitino dark matter is permitted. 

Turning now to panel (b) of Fig. d where the choice m3/2 = 100 GeV is made, we see 
a near- vertical black line at mi/2 ~ 250 GeV: the gravitino is the LSP only to its right. 
The tnsp = 10^ s line has disappeared to larger mi/2, and is not shown. In this case the 
^3/2^^ constraint is much more important than in panel (a), forcing mo to be relatively 
small, simply because m3/2 is larger. The only region allowed by the light-element constraint 
on NSP decays is in the bottom right-hand corner, in the region where the fi is the NSP. 

In panel (c) of Fig. ^ for m3/2 = 0.2mo, there is also a black line to whose right the G 
is the LSP, which is now diagonal, and the VL^i2h'^ constraint is similar to that in panel (b). 
Most of the region allowed by the light-element constraint on NSP decays is in the region 
where the fi is the NSP, though a sliver of parameter space runs above the dotted curve. 

Finally, in panel (d) of Fig. ^ where now m3/2 = mo, the G constraint is more powerful, 
as is the ^3/2/^^ constraint, and the region finally allowed by the light-element constraint on 
NSP decays is again in the fi region. 

Fig. 121 displays a similar array of (mi/2,'n^o) planes for the case tan/? = 35 and /i < 0. In 
the case where m^^ = 10 GeV, shown in panel (a), the most significant change compared with 
panel (a) of Fig. Q is that the 6 — 57 constraint is more important, whilst the VL^i2h^., NSP 
''This line would disappear to larger TO1/2 already for — 20 GeV. 
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Figure 1: The (mi/2,mo) planes for tan/3 = 10,/i > and the choices (a) m-i/2 = 10 GeV, 
(h) m^i2 = 100 GeV, (c) = 0.2mo and (d) = rriQ. In each panel, we show 

nih = 114 GeV calculated using FeynHiggs [24], as a near-vertical (red) dot-dashed line, 
the region excluded by b —>■ s'j is darkly shaded (green), and the region where the NSP 
density before decay lies in the range 0.094 < fi^^p/i^ < 0.129 is medium shaded (grey- 
blue). The (purple) dashed line is the contour where gravitinos produced in NSP decay have 
^3/2^^ = 0.129, and the grey (khaki) solid line (r = 1) is the constraint on NSP decays 
provided by Big-Bang nucleosynthesis and GMB observations. The light (yellow) shaded 
region is allowed by all the constraints. The contour where = nif^ is shown as a (red) 
diagonal dotted line. Panels (a) and (c) show as a black solid line the contour beyond which 
tnsp < 10'^ s, the case not considered here. Panels (b), (c), and (d) show black lines to 
whose left the gravitino is no longer the LSP. 
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decay and tnsp constraints do not change so much. The net result is to leave disconnected 
parts of both the x ^^nd fi regions that are allowed by all the constraints. 

The most obvious new feature in panel (b) of Fig.|2lis the rapid-annihilation funnel, which 
affects both the ^3/2/^^ and NSP decay constraints. The former acquires a strip extending to 
large mi/2 and mo, whereas the latter would have allowed a region at large mi/2 ^ 1500 GeV 
that is excluded by ^3/2/^^- Combining this and the NSP decay constraint, we again find 
two disconnected allowed regions, one in the x NSP region and one that is almost entirely 
in the fi NSP region. 

The rapid-annihilation funnel is also very apparent in panel (c) of Fig. |2l which displays 
the case m^^ = 0.2mo, where again a strip allowed by both the ^3/2/1^ and NSP decay 
constraints extends to large mi/2 and rriQ. There are again disconnected allowed regions in 
the X and (mainly) the fi NSP region. Note that this is constrained at large mi/2 and small 
mo by the t^sp constraint. Finally, in panel (d) of Fig. |2l for m^^ = mo, the region allowed 
by the G LSP, VL^/2h^ and NSP decay constraints is restricted to the part of the (mi/2,mo) 
plane where the fi is the NSP. 

Fig. El displays a similar array of (mi/2,mo) planes for the case tan/3 = 50 and /i > 0. 
The general features of the planes have some similarities to those for tan/3 = 35 and < 0. 
There are differences in the interplays between the (^3/2/^^ and NSP decay constraints, but 
an important difference is the relative weakness of the 6 ^ 57 constraint. This has the 
consequence that allowed x and fi regions are connected for tan/? = 50 and /i > 0. It is 
interesting to note that this is the only case where the putative constraint imposed by the 
muon anomalous magnetic moment impinges on the allowed region, as shown in panels 
(a) and (c). 

We have seen in the above examples that many of the allowed parts of the (mi/2,mo) 
planes are confined to regions where the NSP is the fi. 

6 Conclusions 

We have analyzed in this paper the possibility of gravitino cold dark matter within the 
CMSSM framework. Combining accelerator and cosmological constraints, particularly those 
from h s'-f, VL^i2h? and the light-element constraint on NSP decays, we have found allowed 
regions in the (mi/2,'n^o) planes for representative values of tan/? and the sign of /i and 
different values of m3/2. Standard calculations of the NSP density before decay based on 
freeze-out from equilibrium yield allowed regions where either the lightest neutralino x o^' 
the lighter stau fi may be the NSP. 
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= 10 GeV , tan p = 35 , n < mj/j = 100 GeV , tan |3 = 35 , n < 




mi/2 (GeV) mi/2 (GeV) 

Figure 2: As in Fig. for tan/3 = 35 and ^ < and the choices (a) 7713/2 = 10 GeV, (b) 
"^3/2 = 100 GeV, (c) m^j2 = 0.2mo and (d) 7/13/2 = mo. The light (yellow) shaded regions 
are allowed by all the constraints. 
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m3/2 = lOGeV , tan (3 = 50 , n > m3/2 = lOOGeV , tan |3 = 50 , n > 




mi/2 (GeV) mi/2 (GeV) 

Figure 3: As in Fzg's. Q an(i0 for tan f3 = 50 and fi > and the choices (a) 1713/2 = 10 GeV, 
(b) m3/2 = 100 GeV, (c) 7713/2 = 0.2mo and (d) 7/13/2 = mo- In addition to the quantities 
plotted in the earlier figures, here we also plot grey solid lines where = 44.5 x 10^^^, which 
cut off at small rriQ the allowed regions in panels (a) and (c). The light (yellow) shaded 
regions are allowed by all the constraints. 
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One limitation of our analysis is that it is restricted to tnsp > 10^ s, in order to avoid 
issues related to the hadronic interactions of NSP decay products before they decay. Also, 
in this paper we have not discussed at much length what part of parameter space may be 
allowed in the focus-point region. Finally, we have analyzed here only a few examples of the 
possible relationship between 7713/2 and the CMSSM parameters ttlq and mi/2- 

For these and other reasons, there are still many important issues to analyze concerning 
the possibility of gravitino dark matter. We have shown in this paper that such a possibility 
certainly exists, and that the allowed domains of parameter space are not very exceptional. 
We consider that gravitino dark matter deserves more attention than it has often received 
in the past. In particular, this possibility should be borne in mind when considering the 
prospects for collider experiments, since the allowed regions of the {mi/2,mo) are typically 
rather different from those normally analyzed in the CMSSM. Vive la difference! 
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